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A theory of fluctuations in two-band superconductor MgB2 is developed. Since the standard 
Ginzburg-Landau (GL) approach fails in description of its properties, we generalize it basing on the 
microscopic theory of a two-band superconductor. Calculating the microscopic fluctuation propa- 
gator, we build up the nonlocal two-band GL functional and the corresponding time-dependent GL 
equations. This allows us to calculate the main fluctuation observables such as fluctuation speciflc 
heat and conductivity. 
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I. INTRODUCTION 

The celebrated phenomenological Ginzburg-Landau 
(GL) theory of superconductivity was introduced in the 
seminal paper by Ginzburg and Landaui and has been 
proving ever since to be one of most fruitful and uni- 
versal tools in the description basic properties of super- 
conductors. The examples of its success range from the 
prediction of the Abrikosov vortex stated to the recent 
advance in understanding the complex vortex phase dia- 
gram of the high-temperature superconductors^ and the 
properties of mesoscopic superconducting systems^-. The 
Gor'kov's derivation of the GL equations from the BCS 
theorji^ that have put the GL theory on a firm micro- 
scopic basis and related phenomenological constants to 
material parameters, completed the construction of the 
GL theory of superconductors. 

The GL theory describes well the properties of almost 
all superconductors near transition temperature and is 
successful even when dealing with the superconductors 
with quite complicated band structures. The notorious 
exception from the rule is recently discovered multiband 
superconductor magnesium diboride (MgB2). As it was 
shownSjZ, the GL theory applies to MgB2 only in the 
very immediate vicinity of the transition temperature, 
Tc, i.e. within the interval which turns out to be much 
more narrow than the usual one given by the condition 
|r — Tel ^ Tc- The origin of this narrowing of the in- 
terval of validity lies in the sophisticated band struc- 
ture of the material and reflects the specific interplay 
between its single-electron and superconducting charac- 
teristics. Namely, the Brillouin zone of MgB2 consists of 
the two families of bands: the quasi-two-dimensional cr- 
bands with the strong superconductivity and the weakly 
superconducting three-dimensional 7r-bands. Due to the 
large difference in anisotropy, the c-axis coherence length 
of the TT-bands, £,2,z^ is much larger than the c-axis coher- 
ence length of the cr-bands, S_i,z- Formation of the global 
coherent superconductivity with the unique order param- 
eter implies the appearance of the associated unique ef- 
fective coherence length S,z (T), which turns out to be 
much smaller than ^2,2 at almost any temperature. As 
a consequence, the GL applicability interval shrinks to 



the parametrically narrow range of temperatures where 
S,z (T) ^ S.2.Z- Beyond this temperature range the system 
is strongly nonlocal along the z-axis; to describe such a 
nonlocality one has to employ a generalized nonlocal GL 
model^. One of the spectacular manifestations of the 
non-GL behavior in MgB2 is the strong temperature de- 
pendence of the Hc2 anisotropy close to T^. 

The description of superconducting fluctuations is one 
of the major fields for the application of the GL the- 
ory. Since the standard GL approach turns out to be 
insufficient for MgB2, a generalization of the GL theory 
is required in order to describe its fiuctuation proper- 
ties. The variation of the order parameter on the scales 
smaller than the largest intrinsic coherence length means 
that the usually assumed local approximation does not 
hold anymore and the corresponding short-wavelength 
fluctuations have to be taken into account. 

In the present paper we develop a nonlocal theory of 
superconducting fluctuations that applies to a strongly 
anisotropic two-band superconductor^ and show that 
the short-wavelength fluctuations are essential in it near 
the critical temperature. We start with the derivation 
of the microscopic fluctuation propagator for a such two- 
band model. Then we use it to nonlocal two-band GL 
functional and the corresponding time-dependent GL 
(TDGL) equations. This, in particular, allows us to cal- 
culate the kinetic and thermodynamic observable quanti- 
ties including the fluctuational specific heat and conduc- 
tivity. 

As we have already stated, the main source of the non- 
GL behavior is the nonlocality in the c direction, i.e., 
the strong wave-vector dependence of c-axis coherence 
length. The conventional local form of the GL equa- 
tions turns out to be valid only within the narrow inter- 
val of temperatures, |T - Td/Tc < Ci,z/'?2,z + ^12 < 1, 
where 512 <C 1 is the relative interband interaction con- 
stant which will be specified below. Beyond this interval, 
the superconducting correlations in the 7r-band become 
nonlocal and their contribution to the effective coherence 
length rapidly decreases. Far away from Tc the effective 
c-axis coherence length is determined only by the cr-band. 
In other words, the effective c-axis coherence diverges for 
T — > Tc faster than it could be expected from the naive 
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GL extrapolation, (,z{T) cx — Tc, started from high 

temperatures. This also leads to the decrease of the effec- 
tive anisotropy factor T{T) ~ ^x{T) / S,z{T) ■ As a conse- 
quence, the temperatures dependencies of all fluctuation 
corrections exhibit the characteristic crossovers between 
the dominating cr-band regime (far away from Tc) and 
the "true" GL regime (very close to Tc). For example, 
the c-axis component of the paraconductivity diverges 
faster in the immediate vicinity of Tc than one could 
expect from high-temperature extrapolation using the 
Aslamazov-Larkin formula^"'^ while the fluctuation spe- 
cific heat and ah component of the paraconductivity di- 
verge slower than the corresponding extrapolations. We 
will obtain the temperature dependencies of these fluc- 
tuation corrections. 



II. CRITICAL TEMPERATURE AND 
FLUCTUATION PROPAGATOR 



A. Cooper pairing in two-band model 



The BCS theory was generalized to the case of the two- 
band electron spectrum long time agoi^*i& and has been 
extended recently to include the specific features of mag- 
nesium diboride in Refs. ll7lll8llT9ll20ll2lll23 . We briefly 
overview this theory rewriting it in terms of Green's func- 
tion formalism. The two-band BCS Hamiltonian is given 
by 



(p) V'l,p,aV'a, 



p,cr 



E 

.c{.Q.,l3.(7,a' 



(1) 



where p a ■0q,p,o- ai'e the creation and annihila- 
tion fleld operators in the Heisenberg representation for 
quasiparticle in band a with momentum p and spin a, 

C« (P) = (p- PFa) (2) 

is the quasiparticle spectrum, v^, ppa a-re the Fermi ve- 
locity and momentum of the a-band. The matrix na- 
ture of the electron-electron interaction —gap in O 
fleets the possibility of the interband interactions. The 
free electron Green's functions for each band have the 
usual form Gq {r,r',T,T') = —i (Tril'a'4'2): where Tr is 
the time-ordering operator. In the Matsubara represen- 
tation, 

Ga (P, £n) = {i£n " Cq (p))"^ (3) 

with £„ = 2TiT{n + 1/2) being the fcrmionic Matsubara 
frequencies. 

Now we turn to calculation of the fluctuation propaga- 
tor Lap which characterizes the properties of fluctuation 
Gooper pairs and their effect on observable quantities of 
superconductor above Tc»ii Note first of all, that the in- 
terband electron interactions do not result in the Cooper 
pairing of the electrons from different bands but rather 
lead to the transfer of the pairing correlations between 
the bands. Indeed, the Cooper pairing means the ap- 
pearance of superconducting-type correlations between 
two similar states obeying the condition of the time re- 
versal symmetry. This means that pairing is possible 
only for electrons belonging to the same band, otherwise 
the electron states are too diverse (in terms of the plane 
waves description, their momenta are not the opposite) 
and the integral of the product of their wave functions is 
zero. Thus the formation of unique condensate should 




FIG. 1: Graphic representation of Dyson equation Q for 
the fluctuation propagator L^^ (wavy line). The Greek let- 
ters indicate the band indices. The black dot represents the 
coupling-constants matrix —Qap, the loop represents the po- 
larization operator 11^,4, and the shaded triangle represents 
the scattering vertex C-y-y. 

be understood as the result of the intraband electron 
correlations and subsequent "travel" of the Cooper pair 
from one band to another due to off-diagonal interaction 
components gi2 and g2i- Hence, in terms of diagrams, 
the entrance and exit lines of the fluctuation propaga- 
tor must belong to the same bands (see Fig. ^) and it 
can be presented as the 2x2 matrix Lap (q, ^k) where 
ilk = 27rTfc are the bosonic Matsubara frequencies. The 
corresponding Dyson equation for fluctuation propagator 
can be written in the ladder approximation as (see Fig. 

m 

Lap = —gap + gaj^jSLsp, (4) 

where Hap {ci,i~lk) is the matrix polarization operator 
which is determined by two Green's functions loop^i. 
Without interband electron scattering the only nonzero 
components of the operator Lap are the diagonal ones. 

Turning to the role of electron scattering, let us note 
that at the flrst sight the short-range impurity potential 
might lead to the scattering of electrons all over the whole 
Fermi surface which would result in the equalization of 
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all gap values and reduction of Tc- However, it is well es- 
tablished now that this does not happen in MgB2: even 
in the samples with rather strong intraband scattering, 
interband scattering remains rather weak. This property 
gave possibility to fabricate samples with very high up- 
per critical fields (up to 50 tesla) which have transition 
temperatures only slightly smaller than clean materiaP^. 
Weakness of the interband scattering has been explained 
by Mazin et al. in Ref. ^3 who argued that different par- 
ity symmetry of the a- and tt- orbitals leads to strong 



reduction of the impurity scattering matrix element be- 
tween these bands. For simplicity, we will completely 
neglect the interband scattering. In this case the off- 
diagonal components of the polarization operator van- 
ish. The averaging over impurities position will result 
just in the usual renormalization of the Green's func- 
tions: En — >■ £n + l/2TQ (tq is the corresponding intraband 
scattering time) and to the appearance of the scattering 
vertex part Caa in the expression for the polarization 
operator 



Haa (q, f^fc) ^Ty^ i J^^Caa{(l,en+k,-£n)Ga (p + q, £„ + ^k)Ga (^P, Sn) ■ 

^ J {2tt) 



(5) 



The scattering vertex part can be calculated in the ladder 
approximation^. In the case of a dirty metal {Traa ^ 1) 
this givesii 



Gaai.^: £n-\-k : ^n) - 



' (e„+/c£„) , (6) 



leading to the following result for the polarization oper- 
ator 



In — ^= — ib { — 
27rT ^ I 2 



\^k\+Dc.q^ 
AttT 



(7) 



where Vq, is the density of states in the band a, is the 
Debye frequency and ip {^) is the digamma-function. We 
introduced above the notation for the diffusion-coefficient 
tensor acting in the momentum and band spaces: 



(8) 



The diffusivities Da. a determine the band coherence 
lengths as ^ = nDa,a/&T. For magnesium diboride 
the ratio 



(9) 



is large, r ^ 1, due to the large difference between the 
band Fermi velocities in the c directioni2ii&. This param- 
eter will play an important role in the following consid- 
eration. 

The inversion of the Dyson equation Q gives 



n. 



(10) 



We will use now this general expression to reconsider the 
main properties of two-band superconductivity and to 
describe the corresponding fluctuation properties. 
B. Critical temperature 

The critical temperature of a two-band superconductor 
is determined by the condition det = taken at fifc = 
and q = 0: 

/ .922 , '^1E^d\( 511 , '^1E^d\ .921512 „ 

-.iln^;^j( ^-.,ln— ^ 1—^ = 0. 



\dct 5 



■nTc 



ttTc ) det^5 



with ln7£; = Ce ~ 0.577 being the Euler constant. In- 
troducing the coupling-constants matrix as 



V\g2\ J^2ff22 



we can find that the transition temperature is determined 
by its largest eigenvalue, 



A = ^ + y — + A12A21 

with A± = All i A22, and it is given by the BCS-type 
equation: 

ln?2^ = A-i. 
ttTc 

Following Ref. 0, we introduce the inverse coupling- 
constants matrix 
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where / is the unit matrix and det A = A11A22 — A12A21. 
From the definition^ of the matrix W is evident that it 
is degenerate, detVF = W11W22 — W12W21 — 0. Having 
in mind apphcations of our theory to MgB2, we wih use 
the numerically computed effective coupling constants for 
this compound^: 



^6.7 



Wn W12 
W21 W22 



0.088 
0.424 



0.535 
2.56 



(11) 



We see now that apart from the already mentioned large 
ratio r Q, another characteristic small parameter, the 
relative interband coupling. 



^12 



Wn 

W- 



0.035 < 1 



22 



(12) 



appearsiiii. In what follows we will demonstrate that it is 
the interplay between these two parameters that defines 
the rich picture of fluctuations that occurs beyond the 
traditional GL theory. 



C. Fluctuation propagator 



In terms of the notations introduced above, the inverse 
matrix for the propagator HlOfl assumes the form 



Wu+e + f3i{\nk\,q) 

-W21 W22 



-W12 
■e + (32i\nk\,q) 



lyi 

1^2 



(13) 



where e = lii{T/Tc), 



and 



8T' 



f3{x) EE V (1/2 + 2x/7r2) - ^ (1/2) 



X, a; <C 1 

Inx + Cij +ln(8/7r2), a; > 1 ' 



(14) 



Note that due to the identity W12V2 — ^211^1, the matrix L ^(Jl/j, q) is symmetric. 
Calculating detL~^ and inverting the matrix, one can present L(rife,q) as 



W22 + e[e + l32{\nklq)] _ 
W12 Wii 



L(r!fc,q) = ^ 

e + Wii + Oe 

where we introduced the following notations 

Wiu = 



W21 

[6 + /3i(|r!fc|,(7)] 








1^0 



P2{\nklq) + 



22 



■0/32(|Ofc|,g) /3i(|0fc|,g) 



(15) 



Wii + W 



1 



22 



Wii + W22 



Since we are interested in the most singular contribution with respect to e, we omitted the terms and e as compared 

to W22- 

In order to derive microscopically the TDGL equation and study the effect of fluctuations on transport phenomena 
one has to perform the analytical continuation of the propagator (|15|l from the set of the imaginary Matsubara 
frequencies i^k = iniTk, fc = 0, 1, 2 . . . to the upper half-plane of the complex frequencies. Since the function ip^z) 
has poles at points z„ — 0, —1, — 2 . . . the analytical continuation to the real frequencies 17 can be performed by the 
simple substitution flk —ifl in the argument of the hmction (^^.^^'^^ When this argument is small {il, Daq^ <C T) 
the function f3 is reduced to the Fourier transform of the diffusion operator: 



P ^{\^k\+D^q^) 



. 7r51 TiDaq^ 



8T 



8T 



Performing the analytical continuation, using W22 ~ 1, and accounting for the small parameters of the model {Wn 
S12, 0, e, ri/T, Diq^/T ^ 1), we obtain the analytically continued propagator L^{q,il) in the form 



L«(q,f7) 



1 



e + pi {-in,q) + 



W12 



W21 



S'l2+9[e+/3i(-in,g)] 



Si2/32(-»0,g) ^ _^ 

l+e/32(-iO,g) \ i>i(l+e/32(-iO,(;)) 1^2 (l+9/32(-in,g)) 



(16) 
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Below we will analyze the behavior of this matrix in dif- 
ferent temperature intervals. 

1. Ginzburg- Landau regime 

Let us take the reduced temperature so small that both 
/?! and (32 functions can be expanded. At the end of this 
subsection we will arrive at the analytical criterion for 
this condition to be satisfied. In this case the propagator 
significantly simplifies (we use here the definition 



1 / J/l^l W2XIV2 



(17) 

where ^gl = 7r/8T, and the effective coherence length 
components is given by 



C = W^llCl.a + W22ei,a « U + ^1 



^2, a 



(18) 



Typically S^2,x ~ S,i,x and, due to the presence of a small 
coefficient 5*12, the contribution from the 7r-band to the 
transversal component of the effective coherence length 
(a = X, y) turns out to be small and can be ignored. On 
the other hand, the c-axis motion in the 7r-band, due 
to the large value of the ratio r = £,2,z/&,z can play 
an important role for the longitudinal component of the 
effective coherence length 

g = (1 + Si2r) 



and when S12 ^ l/r it may significantly exceed ^. 

The characteristic momentum is determined by the 
diverging GL coherence length: 



(19) 



Therefore the standard GL theory is valid only at tem- 
peratures so close to Tc that z9z ^ 1 ^^^d where both f3 
functions can be expanded. The corresponding condition 
acquires the forr^^ 



el.. 



5*12 = 1/r + S'12. 



(20) 



2. Beyond the GL regime 

Let us return to the propagator (|15(l and analyze its 
behavior in all the vicinity of the transition, e <C 1, per- 
forming allowed expansions and simplifications. The dif- 
fusion in the first band is slow, therefore we can expand 
the corresponding /3-functions. On the other hand, the 



can be ex- 



function /32(-jr2,q) = P fp{-iVL + D2q^) 

panded only with respect to and q\\ and one has to 
keep the full nonlinear dependence on qz- As a result, 
we obtain the expression similar to (|17|l except for the 
appearance of the explicit g^-dependence of the GL coef- 
ficients 7GL 7 {lz) and Q {qz) 



i«(q,f]) = - 



1 



i-f{qz)n + ea {qz)ql 



'^2(1+3/32) 
\ 1^1(1+0/32) 1^2(1+5/32) 



W12 



(21) 



with 



Here 



lilz)^ IGL 



<S'l2 



l + e[32{qz) 



P2iq.) 



S12 



l + 0f32{qz) 



ez {<!.) = ei,z + 



P2{qz) 



i + ep2{qz) ^izi. 



2 ^2,z- 



(i'{x) = dp{x)/dx, 
/32(g.)^/3fe%9.'), 



(22a) 

(22b) 
(22c) 



(23a) 
(23b) 
(23c) 



One can see that when the condition (|20|) is satisfied, 
the characteristic qz (see (|19|l ) is small and the expres- 
sions (|22af) . H22b|l . and (|22c|l reproduce the TDGL coef- 
ficients. In the region 



1/r + S'12 < e < 1 



(24) 



the argument ^2 z^i (3 function (|23bp for the es- 

sential momenta qz becomes large and it can not be ex- 
panded anymore. This means that the gradient expan- 
sion needed for validity of the GL regime fails. In this 
region of temperatures the contribution of the 7r-band 
rapidly decreases and the main role passes to the cr-band. 
Typically, the 7r-band strongly contributes to (qz) and 
gives only small corrections to 7 (qz) and (qz)- 
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Now we have in place all the elements required for the 
microscopic calculations of fluctuation effects in a two- 
band superconductor. However, because of the complex 
band structure and necessity to take into account the 
short wavelength fluctuations the diagrammatic calcula- 
tions present itself a bulky calculus. In what follows we 
establish another route to address fluctuation phenom- 
ena. Namely, we will re-derive the GL functional and 
extend the standard GL scheme to treatment of fluctu- 
ations for the two-band superconductor, and then apply 
this modified GL approach to calculations of specific heat 
and paraconductivity. 



III. NONLOCAL TWO-BAND GL 
FUNCTIONAL AND TDGL EQUATIONS 

A. GL functional and TDGL equations 

Knowing the explicit form of the fluctuation propa- 
gator H13|l . one can write down corresponding GL free- 
energy functional Fgl = Fgl + ^gl- "^^^ complete pro- 
cedure of its microscopic derivation is given in Ref. ^3 
and here we will present only the specific expressions for 
GL coefficients corresponding to the model under consid- 
eration. 

The quadratic in order parameter Aq part of the GL 
functional F^q]^ is expressed in terms of the linearized GL 
Hamiltonian density i?Q^(q) = (q, i7 = 0) : 



= / drA;i?„^(q)A^, 



(25) 




and q — —iW — ^A. Here the matrix 



1^1 (Wii + e) -1^2^12 
~viW2i V2 {W22 + e) 



plays the role of the GL coefficient a. In general, we have 
to keep nonlinear dependence on in Hapi^ci) meaning 
that Hafsi^q) in Eq. (^5)1 is not a simple second-order 
differential operator. Note that the similar generalization 
of the GL functional has been performed by Maki2i in 



order to describe a dirty superconductor in the vicinity 
of the Hc2{T) line. 

The coefficients in the fourth order term of the GL 
functional F^q]^ do not change their form with respect to 
the noninteracting bands case: 



(28) 



with 



5, = 7C(3):/„/(87r2T2) = jy^bcL 
and C(3) « 1.202. 

Now one can write down the linearized TDGL equation 
in the formii 

(q, n) Af3=[~in-fo.0 {qz)+H^f3 (q)] A^ = 
with the matrix of TDGL dynamic coefficients 



lapiqz) 



ST 



1^1 
Jy2p^{qz) 



(29) 



which directly follows from the dynamic part of the fluc- 
tuation propagator. We see that there are two essential 
differences from the standard GL approach: (i) strong 
nonlocality along z-direction (ii) two-component charac- 
ter of the order parameter. In the case where the param- 
eter W22 is large, one can reduce the functional with two 
order parameters to the functional for the single band- 
averaged order parameter^ 



B. Spectral properties of L ^(q, f2) 

The TDGL operator (q, n) ^ H (q) - i^j {q^) is 
the 2x2 matrix defined on the band index space. It is 
convenient to diagonalize it. The eigenvalues L~^(q) = 
£m(q) — i^lmiqz) and the normalized eigenstates ipm.fs 
obey the equation 



[Hap (q) - i^-faP (qz)] ^m,fj = -^m^(q, 



(30) 



where m = 1,2 is the mode index. Superconducting in- 
stability corresponds to the vanishing of the m — \ eigen- 
value at t,q,Vt — 0. 

Calculating the determinant of H30|l and equating it to 
zero one finds 



r-w _ 1^1 (Wll + fal(q,^^)) + ^2(Ty22+/32(g.) + fe2(q,0)) 

-'^1,2 iqj "J — 2 



± 



yi {Wii + /ii(q, ^)) - V2 {W22 + P2{qz) + /i2(q, O)) 



ViV2WiiW2\ 
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with 



h^ici,n) = e + — + ei,aql 



/i2(q,l^) 



a. GL regime - In the case of small fi, e and cor- 
responding to the GL region of temperatures, one can 
find the simplified expressions for the eigenvalues of en- 
ergy: 



^r'(q) 



V2VI- 



{W22 + Wii) 



VIWII + V2W22 

VlWii -f V2W22- 



(31) 
(32) 



The expression (I31|l for L|f^(q) reproduces the GL rela- 
tion \n\ while is not singular. 

b. Beyond GL regime - In generally, the value of 
function P2iQz) may not be small. In the absence of that 
the smallness the more general expansions for eigenvalues 
have to be used: 



Li ^(q) ~ V1V2 



(/i2(q)+/32(g.))T4^ii +hi{ci){W22 +/32(g.)) 



vi /ii(q) + 



viWii -t- V2 {W22 + hiqz)) 
Wii{h2{ci)+(32{qz)y 



{W22 + fi2{qz)) 

Vi Wii {V2 W22 +V1W11'' 



(33) 



V2 {W22+ l32{qz)) + ViWi 



V2{W22+f32{qz)) ■ 

(34) 



We will need further the eigenvector for the singular 
mode, m = 1, in the zero order with respect to the small 
parameters ha 



where 



^lAqz 

■01,2 («z 



a{qz) 



^l-aHqz) 
-a{qz) 



y/lV^2WuW22 
V2 (1^22 +/32(gz))' 



(35) 



Now we are prepared to revise the results of the standard 
fluctuation theory and to generalize them for the case of 
a two-band superconductor. 



IV. FLUCTUATION PROPERTIES OF 
TWO-BAND SUPERCONDUCTOR 

Before going to detail calculations of different fluctua- 
tion properties, it is instructive to estimate the relative 
strength of fluctuations in the available two-band mate- 
rial, MgB2. It is characterized by the magnitude of the 
Ginzburg-Levanyuk parameter: 



Gi = 



4^^Agr(r,)r, 



where $0 — 2.07-10 ^ G cm^ is the flux quantum, Aa, and 
are the in-plane London penetration depth and coher- 
ence length, and r(T'c) is the anisotropy parameter in the 
limit T ^ Tc- Making use of the values of parameters 
typical for the clean MgB2 crystals: penetration depth 
Xx = 10~^cm (Ref. li^, coherence length = 10~^cm 
(Ref. 0) and the anisotropy coefhcient T{Tc) = 2.5 (Ref. 
M), we obtain Gi « 1.5 • 10~^, what means that fluctua- 
tions in this compound are weak. This conclusion is not 
so surprising, since MgB2 is known to be a good metal 
with large concentration of charge carriers. Therefore, 
identifying experimentally the contribution of fluctua- 
tions in the clean MgB2 crystals is a challenging task. On 
the other hand, the amplitude of fluctuation is expected 
to be much higher in disordered films or in crystals with 
large number of substitution impurities. We stress, how- 
ever, that the effects discussed in this paper hold only 
as long as scattering does not mix bands and this lim- 
its an applicability of our theory for strongly disordered 
materials. Another complication is that increasing dis- 
order in magnesium diboride is usually accompanied by 
doping of the cr-band leading to modification of material 
parameters (e.g., decreasing the cr-band anisotropy)^. 



A. Specific lieat 

We start with the calculation of the fluctuation contri- 
bution to the free energy of the two-band superconductor 
above the critical temperature. It is determined by the 
partition function Z: F = —TlnZ, which, in its turn, 
can be expressed in terms of the determinant of the GL 
matrix Hamiltonian l|26|l 



F = -rin DA DA*exp 



1 



TJ (2.) 



{2T:y det H^p 



where A is an insignificant dimensional constant. Sepa- 
rating the most singular fluctuation contribution (see Eq. 
(|21|l ). one can find 



^™^--T/-^ln 



(2 



(l + 0p2{qz)) {e+ei,aQl + Si2P2{qz)) 



V 



with P2iqz) = P2{qz)/ (^1 + OI32{qz)j ■ The corresponding 
contribution to the specific heat is given by 



T d'^F 
V dT^ 



d\ 



1 



^^ + eiWa + Si2Uqzf 



CO 



dqz 



^+ei,zqi+si2i32{qzr 



(36) 
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Note that the denominator of the the logarithm argu- 
ment coincides with the denominator of the fluctua- 
tion propagator 1)21(1 . Introducing the reduced variable 
u — ^i.zQz/y/^, we rewrite this result in the form conve- 
nient for numerical evaluation: 



C' = 



(37) 



It weakly depends on temperature and has the following 
asymptotics 



where the dimensionless function K(e) = K{e,r, Si2,0) is 
defined as 



2du 



TT l+u2 + (5i2/e)/32(r-eu2) 



(38) 



^'^^ =^i,z/^z, fore«l/r + 5i2 

{Si2/e) In [Cr (e - 5i2)], for 1/r + ^12 « e « Su HrSi2) 
1 - for S12 In(r5i2) « e « 1 



(l + 5i2r) 



2£ l+ein(Cire) ' 

with Ci = Wtt-'^je exp (-2) « 0.391. 

The formulas ((37|I - H38|I work in the entire region e 1. The intermediate asymptotic in k appears only in the case 
rSi2 ^ 1 which is valid for MgB2»ii One can see that the main difference between the GL and non-GL regions in 
the temperature dependence of the fluctuation heat capacity correction is the change in the coefficient of the e~^/2 
dependence in factor S,z/£,i.z- 



B. Paraconductivity 



The paraconductivity in the phenomenological GL approach can be expressed via the eigenvalues of the GL Hamil- 
tonian and the matrix elements of the "velocity" operator = dHaa (q) /dqa- The approach of Ref. can be 
extended it to the case of two-band model. The general formula for the paraconductivity tensor then is: 



h 

4Te2 



(27r)^ ^ £n£k {ln£k 



■ Snlk) 



2el 



72^2^^22 

2sl 



r 



2^?2^^21 



£i£2 (71 ^ei + 72 ^£2) 



(39) 



with v^f, = J2a '^n,aV%i^k,a (we rcstorcd dimensional 
units in this formula). Let us stress that the summa- 
tion here is performed over the mode indices rather than 
the band indices in other parts of the paper. 

The main contribution to the paraconductivity comes 
from the projection to the singular mode (the first term 
in the square brackets in Eq. H39|l 'l. Keeping only these 
terms and using results H33() for the eigenvalue Lj^^(q) = 
ei(q) — iVt'^i{qz) and ((35|l for eigenvector, we derive 



ire 



2 r d^q ^l,all '^ + Sl2raP'a 



^'^^^^ fe + ^2 2^512/32'' 



(40) 



/3y (1 + ^/32 andr, 



with i3[ = 1/(1 + ^/32)', P'2 - 

Cla/Cia = Below we evaluate the in-plane and 
z-axis components of paraconductivity separately. 

1. In-plane component 

In the case of an in-plane paraconductivity we can ne- 
glect the small renormalization of the in-plane velocity 
in formula H40|l . perform the qn integration and get 



32S7r 



dqz 



1 



e + ^L«.' + ^12/32(9.) 
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1 \ GL region 




V x' intermediate 






o-band 







0.04^0.06 



FIG. 2: Temperature dependencies of the paraconductivity 
components (solid lines). Left panel shows the in-plane com- 
ponent in units of e^/(32ft^z,i) (function K{t)/y/e, see Eqs. 
(I41I I and 1381 1'). This curve also gives the fluctuation specific 
heat in units l/(87r^i_j;^i,z), see Eq. (1371 . Right panel shows 
the z-axis component in units of e^^z,i/(32fi,5^_i) calculated 
from Eq. 1421 . Both components were calculated using pa- 
rameters = 0.035, 9 = 0.377 and r = 300. For comparison 
we also show the Ginzburg-Landau asymptotics (dotted lines) 
and (T-band contribution to paraconductivity (dashed lines). 
One can see that in-plane component diverges slower and z- 
cods component diverges faster than the paraconductivity in 
a single-band anisotropic superconductor. In the a"z-plot we 
indicate three characteristic regions with different functional 
dependencies of fluctuation parameters. 



One can see that the in-plane conductivity has exactly 
the same temperature dependence as the fluctuation spe- 
cific heat and therefore can be represented in the form 
analogous Eq. (|36() : 



= K(e), 



(41) 



where the function K(e) is defined by Eq. (|38l) . The 
above results show that a'^ has a 3D character but di- 
verges slower than the specific l/\/e-law. Numerically 
calculated dependence a'^ (e) with parameters typical for 
MgB2, = 0.035, B = 0.377 and r = 300, is shown 
in the left panel of Fig. [21 together with the GL asymp- 
totic for T — > Tc and the single tr-band curve. Taking 
the estimate for the coherence length of MgB2 crystals 
asS ^1,2 ~ 2nm, we find the typical scale for the fluc- 



e7(32Sa,^), as ct: 



40 [n 



^. This scale must be much higher in dirty MgB2 
The in-plane paraconductivity in MgB2 films has 



tuation correction 
cm]~ 
films 

been studied recently in Ref. but in the analysis of 
the experimental data the two-band nature of MgB2 was 
not taken into account. Nevertheless, it was found that 
the paraconductivity indeed diverges slower than 1 / y^, 
in agreement with our results. 



2. z-component 

The longitudinal component of paraconductivity a'^ is 
given by the Eq. (|40|) . Here the velocity renormaliza- 
tion turns out to be essential and it cannot be omitted. 
Performing integration with respect to and using the 



introduced above reduced variable u, the result can be 
rewritten as 



4du 



(42) 

In the Ginzburg-Landau regime, e <C l/?'+5'i2, we obtain 

, ^ e^a.^yi + '5'i2r 
32^2,7? ■ 

In the regime of nonlocal fluctuations (e 1/r + 5i2), 
out of the GL region, the behavior cr^ becomes quite pe- 
culiar. Essential contributions to the integral in Eq. 14U|) 
occur from two regions of qz'. from ^ 1/^2, z in ar- 
guments of functions P2 and (3'2 (mainly 7r-band contri- 
bution) and from ~ V^/S.i,z (mainly cr-band contri- 
bution). This corresponds to the ranges u 1/re <C 1 
and u ~ 1 in the reduced integral Eq. H42|l . We evaluate 
separately contributions from these ranges, or, what is 
the same, from tt- and a- bands in the most interesting 
case rSi2 ^ 1 relevant for MgB2. 

In the range u ^ 1/re one can drop cr-band terms in Eq. 
(I42|l . As we consider the regime e ^ S'12 we can neglect 
the term with /3 in denominator. Dimensionalizing the 
integral, one can obtain the following result for this term 



e'Ci,zS'i 2V^ 
32ftCL 



C2i0), 



4dv 



gives C2(0) 



(43) 



1.185 and 



Numerical calculation 
C2(0.377) « 0.235. 

In the range u ^ I the main contribution occurs from 
the cr-band and the terms proportional to S12 (7r-band 
terms) can be treated as small perturbations. Expansion 
with respect to these terms leads to the following result 
for the paraconductivity contribution 



with 



16 

TT 



du- 



S12 
~2^' 



(44) 



(1+^2)^ 

In (Care) 



er (32 



/32 



with C3 = IQn-'^-fE exp(-4) sa 0.053 and C4 « 2.9. Let 
us note that in this range of temperatures the contribu- 
tion from the 7r-band comes with the negative sign, as in 
the cases of the heat capacity and in-plane paraconduc- 
tivity. This means that the 7r-band contribution changes 
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sign with increasing temperature. Comparing contribu- 
tions (|43|l and (|44l) . we find that the 7r-band term domi- 

1/3 

nates in the interval of temperatures e < (^'S'12) 



Therefore in the case rS'12 ^ 1 the z axis paracon- 
ductivity cr^ has three asymptotic regimes: 



(l + 5i2r)i/% fore«l/r + 5i2 
■^9!^ W^^i^C^^W/e^/^, forl/r + 5i2«e« f^:^ . (45) 

I for max[l/r + 5i2,5i2(r5i2)^/^]«6«l 



In the case rS'12 < 1 the intermediate asymptotic disap- 
pears. 

The right panel of Fig. |2 shows the numerically cal- 
culated dependence of a'^ (e) with parameters specified 
in the captions. For comparison, we also show the GL 
asymptotics and the cr-band contribution. 



V. FINAL REMARKS 

In this paper we have presented a microscopic deriva- 
tion generalizing the conventional GL description onto 
two-band superconductors. We have further applied the 
developed approach to the investigation of the fluctu- 
ation phenomena in MgB2. The important feature of 
our approach is that the derived nonlocal GL functional 
takes into account not only the long-wavelength fluctua- 
tions (as is the case of the conventional GL theory), but 
also the short-wavelength fluctuations, which enables to 
significantly extend the range of validity of the GL tech- 
nique. This approach not only permits the study of the 
thermodynamic characteristics of the superconductors, 
but provides a convenient technique for calculation of 
transport coefficients. 

It is necessary to underline that we succeeded to 
solve analytically the problem of accounting for short- 
wavelength fluctuations due to the specific for our model 
large difference in the intra-band diffusion coefficients 
Da- Fortunately, this assumption corresponds to the 
practically interesting case of magnesium diboride. 

One can see that the main qualitative effect that de- 
termines the unusual behavior of fluctuations in the 
anisotropic two-band model consists in globalization of 
superconductivity in the immediate vicinity of transition 
temperature (e <C 1/r + Su). It occurs due to the possi- 
bility for electrons of both bands to participate in fluctua- 
tion pairing and to exchange of fluctuation pairs between 
them. Formally, this process manifests itself in the ap- 
pearance of long range superconducting correlations on 
the scale ^ £.i,z- 

Turning the role of fluctuations in MgB2, we refer 
to the formulas for paraconductivity in a one-band 3D 



anisotropic superconductor—: 

/ ^ I 6 ( AP,^ 

We see that the growth of the effective £,z results in the 
noticeable increase of the z-axis paraconductivity in the 
immediate vicinity of critical temperature (e ^ 1/ r-l-S'12) 
with respect to its high temperature (e > Sf^'^r^^^) ex- 
trapolation based only on the cr-band ffuctuation pairing. 
The crossover between these two 3D regimes occurs in the 
narrow interval of temperatures l/r-|-5'i2 ^ e > Si2'^r^^^ 
where the boson degrees of freedom corresponding to the 
pairings in 7r-band rapidly freeze out, leading to the fast 
decrease (~ e~^) of z-axis paraconductivity. We want to 
stress that this temperature dependence appears due to 
the noticeable contribution of short-wavelength fluctua- 
tions in this range of temperatures. 

In the case of the in-plane component of paraconduc- 
tivity and the fluctuation part of heat capacity, the sit- 
uation is the opposite. The coherence length appears 
in the denominator of the corresponding expressions, so 
tr^ and C" turn out to be suppressed in the vicinity of 
transition with respect to their extrapolation formulas 
from the high-temperature behavior. Similar situation 
occurs in the temperature dependencies of the in-plane 
upper critical field and of anisotropy of the upper critical 
field&Z. 

The obtained results coincide with that ones of the di- 
agrammatic approach but the proposed GL description 
has an advantage of being more physically transparent, 
economic, and universal. The fact that the results de- 
rived via the GL equations coincide with those of our 
microscopic consideration is a convincing cross-check of 
the phenomenological description. The same approach 
can be used to derive other fluctuation properties such 
as the magnetic susceptibility and the field dependen- 
cies of conductivity and magnetization in the vicinity of 
transition. 
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